Probing quasiparticle dynamics in Bi2Sr 2 CaCu 2 08+<5 with a driven Josephson vortex 
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We show that the flux-flow transport of the Josephson vortex lattice (JVL) in layered high- 
temperature superconductors provides a convenient probe for both components of quasiparticle 
conductivity, a c and a a b- We found that the JVL flux-flow resistivity, pff, in a wide range of 
magnetic fields is mainly determined by the in-plane dissipation. In the dense lattice regime (B > 1 
T) pff(B) dependence is well fitted by the theoretical formula for that limit. That allows us to 
independently extract from the experimental data the values of a c and of the ratio a a b/ (o'c'y 4 ) ■ The 
extracted temperature dependence cr a j,(T) is consistent with microwave data. The shape of the 
current-voltage characteristics is also sensitive to the frequency dependence of a a b and that allows 
us to estimate the quasiparticle relaxation time and relate it to the impurity bandwidth using data 
obtained for the same crystal. 



I. INTRODUCTION 

The properties of quasiparticles (QPs) in the high- 
temperature cuprate superconductors are unusual due 
to d-wave gapless pairing in these systems. The con- 
centration of QPs in clean d-wave superconductors van- 
ishes in the limit T — ► because their density of states 
(DOS) has asymptotics p(e) oc e at low energies, e — > 0. 
Such d-wave features of the QPs spectrum has been 
clearly demonstrated by angular-resolved photoemission 
(ARPES) studies (see, e. g., Refs. Q| and confirmed by 
numerous transport measurements. However, QP trans- 
port properties are still in the focus of hot debate and 
new experimental techniques are very valuable in under- 
standing of these properties. In this paper we use the 
driven Josephson vortex lattice to probe both in-plane 
and inter-plane QP transport and compare results with 
other techniques. 

We start with brief overview of the QP physics in the 
high-temperature superconductors. The main issues are 
how impurities modify the QPs spectrum and how they 
affect mobility of the QPs. Two competing effects de- 
termine QP transport at low temperatures. It was rec- 
ognized very early that impurities in unconventional su- 
perconductors destroy the superconducting order param- 
eter in their vicinity increasing the DOS at low energies. 2 
This effect enhances the low-temperature QP transport. 
On the other hand, increase the QPs scattering rate by 
impurities suppresses the QP transport. In addition, 
there is a tendency for QPs localization due to the two- 
dimensionality of superconducting Cu02 layers. 

The low-energy asymptotics of the QP DOS in inho- 
mogeneous d-wave superconductors is a challenging the- 
oretical problem, which is not settled yet. Several ap- 
proximations have been used to calculate p(e) leading to 



different results. In the self-consistent T-matrix approxi- 
mation (SCTMA)ii^ the DOS approaches a finite value 
p(e) ~ p(0) for e < 7, and p(e) oc e for e > 7,, where 
the impurity bandwidth 7^ depends on the impurity con- 
centration and the impurity potential. In the limit of 
strong impurity potential (unitary limit) one gets for the 
impurity bandwidth^ 73 « (fa^Ao) 1 / 2 , where vq is the 
normal-state relaxation rate and Ao is the magnitude of 
the superconducting gap. 

With p(e) known, the QPs transport can be calculated 
in the framework of the Fermi liquid theory. Along this 
line Lee^ predicted a universal (impurity independent) 
low-temperature limit for the intralayer electrical con- 
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(06) 
00 



2 /(7TTO bAo), arguing that both, the 



DOS and the scattering rate are proportional to the im- 
purity concentration and thus cancel in the Drude expres- 
sion for the conductivity. Here n is the carrier concen- 
tration and m a h is the QP in-plane effective mass. In the 
same way universal thermal conductivity was predicted 
by Sun and Maki£ and by Graf et al£ Later it was shown 
in Refs. |3 that electron interactions inside the layers lead 
to the Fermi-liquid corrections to the universal electrical 
conductivity. Durst and Lee 9 found that such corrections 
as well as the asymmetric scattering (i.e. the difference 
between the QPs relaxation rate and the transport scat- 
tering rate) result in the expression for the intralayer elec- 
trical conductivity <J a h(T) = Cqq 6 '*/?. Experimentally ob- 
served values of low-temperature a a b usually correspond 
to (3 > 1, see Ref. ^3 Durst and Lee found also that 
such corrections to the universal thermal conductivity at 
e < 7j are practically absent. 

For the interlayer tunneling conductivity, er c , in highly 
anisotropic layered cuprat.es like BSCCO the authors of 
Ref. ^ have argued that it is universal in the limit of 
low impurity concentration, when electrons tunnel be- 
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tween layers conserving their in-plane quasi-momentum 
(coherently). In contrast to a a t,, a c depends only on 
the intralayer DOS (i.e., on the QPs relaxation rate). 
It is not sensitive to anisotropy of the in-plane scatter- 
ing and in-plane vertex corrections are not important. 

(c) 

In this case the low-temperature conductivity a 0Q is re- 
lated to the Josephson critical current density J c (0) by a 
simple relation similar to the well-known Ambegaokar- 
Baratoff relation in conventional Josephson junctions, 

Coo* = 2esJ c (0)/(7rAo), where s is the interlayer spac- 
ing. In the framework of SCTMA-Fermi-liquid approach 
the temperature corrections to the universal conductivity 
<Jqq at T < 7; were founcUi to be 

v c {T)*<r$[l + {itT)*/18tf\. (1) 

A similar result has been obtained for the in-plane QP 
conductivity neglecting anisotropic scattering and vertex 
corrections £ However, it is known that such corrections 
substantially modify a a b(T). For the thermal conductiv- 
ity, k(T)/T, in the unitary limit, the thermal corrections 
are also quadratic in T/'fi, as for cr c , and only in the 
limited temperature interval linear in T corrections were 
found if the impurity potential is not very strong. 12 

The more elaborate approaches, which take into ac- 
count interference effects, suggest that the QP DOS van- 
ishes at e — > (see recent review in Ref . Il3j) . In particu- 
lar, recent numerical analysis^ have demonstrated that 
QP DOS behaves at very low energies as p(e) oc e a with a 
nonuniversal exponent a. This exponent depends on the 
details of disorder model and particle-hole symmetry in 
the normal state. In the realistic case (binary alloy model 
without particle- hole symmetry) they found a > 0, i.e., 
a DOS suppression at low energies. The energy scale for 
this suppression is given is the resonance energy for an 
isolated impurity Qq with Qq — > as the impurity poten- 
tial increases. This means that at very low temperatures, 
T <C fioi QP transport should vanish. 

Intralayer conductivity was studied by a microwave 
techniquei2*iii*i§iii and infrared spectroscopy>i2ii2i2£ 
These measurements have shown that a a b(T) is not uni- 
versal in the low temperature limit and that its tem- 
perature dependence at T ~ 7.; deviates strongly from 
SCTMA predictions^ see below. The impurity band- 
width can be estimated from the frequency dependence 
of the intralayer conductivity, <7 a b(w). According to 
the model calculation of Ref. 0j the typical relaxation 
rate, 1/r, in a a b(uj) has to be of the order of 7,. Re- 
cent terahertz spectroscopy measurements of a a b{uj) in 
Bi 2 Sr 2 GaCu 2 8+5 (BSCCO) by Corson et alM showed 
that at low temperatures it has a Drude frequency de- 
pendence with the typical relaxation rate 1/t w 1 THz. 
This gives an estimate 7, = 30 — 50K. 

Experimentally a universal value of kqq/T was con- 
firmed by measurements2i*22iS^ of the thermal conductiv- 
ity in BSCCO and YBCO crystals. Particularly, Naka- 
mae et td*L measured thermal conductivity in pristine 
and irradiated BSCCO crystals at T < 0.1 K. Direct com- 



parison of k(T) in pristine and irradiated crystals showed 
almost the same value kqo/T obtained by extrapolation 
to T = 0. However, strong temperature corrections, lin- 
ear in T, to this universal value were found below 0.25 K. 
Such corrections are not anticipated well below 73 sa 20 K 
in the SCTMA approach. One can conclude from these 
measurements, that the upper limit for the low tempera- 
ture behavior, lies below 0.1 K and that temperature 
behavior of k(T) observed so far is in contradiction with 
the theoretical SCTMA predictions. 

Measurements of interlayer conductivity seems to fol- 
low closely to the theoretical predictions of the the 
SCTMA-Fermi-liquid model assuming coherent inter- 
layer tunnelingiii For BSCCO crystals from different 
groups, similar values, « 2 [kfi- cm] -1 were observed 
by measurements of the I-V characteristics, though uni- 
versality of this value was not checked by comparison of 
crystals before and after irradiation as it was done for the 
thermal conductivity^ The temperature dependence of 
<J C (T) at T < ji was found to be in agreement with the 
theoretical prediction, Eq. QJ, and the values of 7$ in the 
interval 24-29 K were derived from fitting, in agreement 
with the data for the relaxation rate r mentioned above. 

The intralayer and interlayer components of the QP 
conductivity are found to have qualitatively different 
temperature dependencies in the superconducting state. 
The interlayer conductivity monotonically decreases with 
temperature decrease in the whole temperature range 
from T c down to lowest temperatures. In contrast, the 
intralayer conductivity has manifestly nonmonotonic be- 
havior: it rapidly increases with decrease of temperature 
in some region below T c reaches a peak at some interme- 
diate temperature, and decreases with further decrease 
of temperature>i2ii£*i£ii£ This behavior is also consistent 
with the heat transport measurements for the electronic 
part of the thermal conductivity^ The peak in the tem- 
perature dependence of a a b appears due to an interplay 
between the temperature dependencies of the relaxation 
rataii and concentration of quasiparticles. The unusual 
increase of the QP conductivity and thermal conductiv- 
ity just below T c is attributed to the fast drop of the 
relaxation rate due to the reduction of the phase space 
for scattering^ Decrease of the conductivity at low tem- 
peratures, below the peak, is caused by a drop in the 
concentration of the thermally activated nodal quasipar- 
ticles. 

Hence, the behavior of QPs in the cuprate d-wave su- 
perconductors is not fully understood yet and new meth- 
ods to probe the QP transport would be useful to resolve 
controversy and provide additional information on the 
characteristic parameters, such as f2o- 

Recently a new method of probing the QP conductiv- 
ities a c and a a b has been suggested^ The method is 
based on the measuring of the losses, associated with a 
transport of the Josephson vortex lattice (JVL)2^ driven 
by steady current across the layers in crystals with in- 
trinsic Josephson interlayer coupling^ The unexpected 
dependence of the interlayer transport on the intralayer 
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FIG. 1: Dense Josephson vortex lattice with period / = 2s in 
the c-direction. Arrows show the direction of currents. 



quasiparticle conductivity in the flux-flow regime is re- 
lated to the spatially inhomogeneous structure of moving 
JVL. Figure n shows the space structure of a stationary 
high-field Josephson latticed Inside the layers the su- 
percurrent oscillates along the direction perpendicular to 
the applied parallel magnetic field (x -axis). The inter- 
layer transport current drives the vortex lattice along the 
x-axes. At small velocities the lattice is practically undis- 
torted. Then the supercurrent inside the layer at a given 
point changes periodically with time, j sx oc sm(2irvt / a) , 
where v is the lattice velocity and a = &o/sH is the 
period of vortex lattice. According to the first London 
equation (see, e.g., Ref. Wftl 



E x = 



c 2 dt 



\ 2 - 

A ab — 



m ab c* 
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an alternating electric field E x with frequency u> — 2?rv /a 
is introduced by a moving vortex lattice. Here n s is 
the density of superconducting electrons. The electric 
field causes relaxation due to quasiparticle current j qx = 
&abE x - Hence, as was shown, both components of the 
QP conductivity a c and a a b contribute to the Josephson 
flux- flow resistivity pft and can be extracted separately 
from the magnetic field dependence of pjf. Moreover, 
the shape of the I-V characteristic related to JVL mo- 
tion at high magnetic fields is sensitive to the frequency 
dependence of the QP conductivity and therefore may be 
used to estimate the QP relaxation time. This approach 
allows one to probe the relatively low frequency range of 
the QP conductivity, 0.01-3 THz, and thus the results can 
be easily compared with the data of microwave measure- 
ments. The goal of the present paper is to demonstrate 
the applicability of this new method for studies of the 
QP conductivity in BSCCO. 

Motion of JVL induced by a steady current across the 
layers results in a specific branch on the I-V characteris- 
tic usually referred to as the Josephson flux-flow (JFF) 
branch. That is characterized by a rapid current increase 
when the voltage approaches a certain limiting value, 
V m , at which the lattice velocity approximately reaches 
the velocity of electromagnetic wave propagation (this 
velocity frequently is refered to as the Swihart velocity 
in analogy with a single long Josephson junction). The 
JFF regime is well known for conventional long Joseph- 
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FIG. 2: (a)Schematic view of junction orientation and ex- 
perimental set up. (b)Normalized temperature dependence 
of resistance along the c-axis of typical long stacked junction 
used in experiment. 



son junctions^ and has also been observed on BSCCO 
mesa structuresi 29 i 30 i 31 Our purpose was to study the 
JFF linear and non-linear regimes on long BSCCO stacks 
at high magnetic fields above 0.5 T when a dense JVL 
is formed. Early experiments on JFF in BSCCO have 
been done mostly at relatively low fields22i2i or at high 
fields of only few fixed values^ In addition, an absence 
of the clear upturn curvature on JFF branch 30 suggests 
significant inhomogeneities in the mesas used in early ex- 
periments. 

We report new measurements that allow us to obtain 
for the first time the parameters 1 jr and 7, for the same 
crystal and correlate the frequency and temperature de- 
pendencies of quasiparticle conductivities. 



II. EXPERIMENTAL 

To obtain high-quality stacks we fabricated samples 
from single crystal whiskers of BSCCO. The thin BSCCO 
whiskers have been characterized as extremely perfect 
crystalline objects^. They grow along the [100] direc- 
tion free of any crucible or substrate and can be entirely 
free of macroscopic defects and dislocations. The stacks 
have been fabricated by the double-sided processing of 
the BSCCO whiskers by the focused ion beam (FIB) tech- 
nique. The stages of fabrications were similar to those 
described in Ref. Figure shows schematically the 
geometry and orientation of the structure with respect 
to the crystallographic axes. We reproduced the overlap 
type of long stack geometry which is known to provide the 
most uniform current distribution along the junction. 28 
The structure sizes were L a — 20 — 30 fim, Lf, = 1 — 2 
jitm, L c — 0.05 — 0.15 /im. Using high-resolution optical 
microscope we selected for the experiment long uniform 
whiskers with a length of 500-1000 fim, a thickness of 0.5- 
1 /im and a width of 10-20 /xm. Four silver contact pads 
have been evaporated and annealed at 450° C in oxygen 
flow before FIB processing to avoid diffusion of Ga-ions 
into the junction body. The fabricated structures have 
been then tested by R C (T) measurements to select ones 



4 



free of inclusions of 2201 or 2223 phases. The presence of 
these phases is indicated as a multiple transition to super- 
conducting state with appropriate drops of R C {T) at 90- 
100K for 2223 phase and at 15-30K for 2201 phase. Only 
single-phase 2212 stacked junctions with a single transi- 
tion at 75-80 K (see Fig. Eb) have been selected for the 
further measurements. A yield of the single-phase stacks 
was quite high, about 30-50%. The oxygen doping level 
of the stacks estimated from p c (T) measurements above 
Tj^ was slightly above optimum, 6 « 0.25. The critical 
current density J c at 4.2 K in the absence of magnetic 
field was 1-2 kA/cm 2 . Measurements of the I-V charac- 
teristics of BSCCO stacks have been carried out in the 
commercial cryostat of Quantum Design PPMS facility. 
The magnetic field has been oriented parallel to the b- 
axis within accuracy 0.1° and has been changed in steps 
of 0.05-0.1 T. In each fixed value of the field the I-V char- 
acteristics have been measured using a fast oscilloscope. 
We have measured 6 samples with similar results. 



III. RESULTS AND DISCUSSION 

Fig. |3 shows the I-V characteristics of a long stack 
#4 (30 x 2 x 0.14 pm 3 ) at T = 20 K at zero magnetic 
field (Fig. G&) and at the field B = 1.5T (Fig. Eb) ori- 
ented along the b-axis. At zero magnetic field the I-V 
characteristics of the stacks show a well-defined critical 
current shown at Fig. [3^ as a vertical trace with the fol- 
lowing switch to the multibranched structure. The I-V 
characteristic is highly hysteretic. The hysteresis loop 
is shown in Fig. |2K schematically by arrows (the first 
switch here corresponds to the jump from the critical 
current to the A th branch). Both features, the hysteresis 
and the multibranched structures of the I-V characteris- 
tics, are well known for Bi-2212 stacked junctions at low 
temperatures i2£ 

In the presence of a parallel field the critical current 
becomes essentially suppressed and the JFF branch de- 
velops. That is characterized by a linear slope at low 
bias, Rff, by a pronounced upturn at higher bias volt- 
ages, and by jumps to the multiple branches at the volt- 
ages exceeding the maximum value V m . The high upturn 
in the nonlinear I-V characteristics in the parallel field 
proves a high quality of our stacks. To get the Rff(B) 
dependence we measured a set of the I-V characteristics 
at some temperatures for the fixed fields increasing step 
by step. At each field we measured Rff- The value of 
Rff has been defined as an extrapolation of the linear 
part of the I-V at \V\ — > 0. As seen from Fig. [3J;, this 
extrapolation can be easily done at high fields when crit- 
ical current is highly suppressed. An accuracy of linear 
extrapolation of the I-V characteristic at | < lOmV 
for definition of Rff was within 5%. Fig. |3t shows an 
evolution of the I-Vs with field. One can see a rapid in- 
crease of Rff and V m with field. The summarized field 
dependencies of Rff and V m for typical sample (#4) are 
shown at Figs. 01 and Before analyzing these data we 



discuss the expected theoretical dependencies. 

Firstly, we will consider the linear limit of the I-V char- 
acteristics corresponding to low JFF velocity. In the sec- 
ond part, we will focus on the high- velocity JFF limit. 



A. Linear flux-flow resistivity. 

The linear flux-flow resistivity of the Josephson vor- 
tex lattice, pff, is determined by the static lattice struc- 
ture and linear quasiparticle dissipation. At high fields, 
B > $ ( 7r 7s 2 ), the Josephson vortices homogeneously fill 
all the layers and the static lattice structure is character- 
ized by oscillating patterns of both c-axis and in-plane 
supercurrents 2 ^ (see Fig.^). Here 7 is the anisotropy ra- 
tio of the London penetration lengths, 7 = A c /A a f,. At 
small velocities this pattern slowly drifts along the direc- 
tion of layers, preserving its static structure. This mo- 
tion produces oscillating c-axis (E z ) and in-plane (E x ) 
electric fields leading to extra dissipation, in addition to 
usual dissipation due to the dc electric field E z . Total 
dissipation per unit volume is given by 



OffEl = a c El 



Et 



Vab I El 



(2) 



where a// = 1/p// is the flux-flow conductivity, (. . .) 
means time and space average, a c — l/p c and a a i, — 
1/Pab are the c-axis and in-plane quasiparticle conduc- 
tivities. An expansion with respect to the Josephson 
current at high fields allows to relate E z and E x with 
E» 

4E Z 

E nz (x,t) = -{-l) n cos (k H x + ujt) , 
IE 

E nx (x,t) = (-1)™— ; 1 sin (k H x + ujt) , 
7/1 

with h = 2TT"fs 2 B/<& and kn = 2-ksB/Qq. Here n is 
the layer index. Finally, we obtain a simple analytical 
formula for the flux-flow resistivrhi 2 ^^ 



Pff(B) 



B 2 



B 2 +B 2 



B n 




v / 27T7 2 s 2 



(3) 



The relative importance of the in-plane and c-axis dissi- 
pation channels is determined by the dimensionless ratio 
r = cr a fc/(cr c 7 2 ). The in-plane dissipation channel dom- 
inates when r ^> 1. 

Experimental dependencies of the JFF resistance, 
Rff(B), on magnetic field B applied along the 6-axis at 
various temperatures are depicted in Fig.0J They are ev- 
idently nonlinear with the type of non-linearity predicted 
by Eq. quadratic growth with the following satura- 
tion. As shown in Fig. 0] the data can be well fitted to 
Eq. © in the wide range of temperatures. Two fitting 
parameters have been used for each curve, i? c , the high 
field saturation value of the resistance, Rff(B), and the 
characteristic field B a , both defined by Eq. Using 
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FIG. 3: The I-V characteristics of stack #4 at 4.2 K without magnetic field (a) and with field B applied along fa-axis, B — 1.5T 
(b), with B increasing from 0.85 T up to 1.5 T (c). The first three branches are not traced on the I-V characteristic in zero 
field and the stack jumps directly to the 4th branch. 



that fit at a set of temperatures we found a temperature 
variation of both parameters. Then, from the tempera- 
ture dependence of R c we can directly get the dependence 
of a c {T), while the temperature dependence B a contains 
information about the ratio a a b(T)/a c (T) for given val- 
ues of 7. For T = 4.2 K we found B a = 3.3 T. We esti- 
mate the value of 7 w 500 at 4.2 K from the value of the 
Josephson critical current density J c (0) = 1.7 kA/cm 2 
and the value of A a fc(0) = 0.2 /im using the well-known 
expressionSi, J c — c$o/(87r 2 s7 2 A 2 h ). That gives a quite 
reasonable value for <r a b at 4.2 K, <r a f,(4.2) = 4- 10 4 (Ohm 
cm) -1 . Finally, we restored the temperature dependence 
of the in-plane quasiparticle conductivity a a b(T) taking 
into account the temperature dependence of 7. We used 
the j(T) dependence extracted from the known data 
for A afc (T)S and A C (T)22 in BSCCO. Actually, 7 2 only 
weakly depends on T, slowly decreasing within 15%, with 
temperature increase from 4.2 to 70 K. 

Figure [fj] demonstrates the temperature dependencies 



of <j a i, and a c extracted from the Josephson flux-now ex- 
periment. As we found, a a b rapidly increases below T c 
with decrease of T, reaches a maximum at about 30K, 
and then drops down at low temperatures. That type of 
Pab{T) behavior has been found earlier in the microwave 
experiments for YBC02& and BSCCQ 15 i 16 i 17 and also 
is consistent with the heat transport measurements of 
the electronic part of the thermal conductivity. 2 - As dis- 
cussed in the Introduction, the peak in the temperature 
dependence of <j a b appears due to an interplay between 
the competing temperature dependencies of the relax- 
ation rate 1 ! and concentration of quasiparticles. We 
found that the value and temperature dependence of 
dab extracted from the Josephson flux-flow experiment 
quite well reproduce the low-frequency microwave data^ 
The sharper increase of a a b{T) below T c may be re- 
lated to weaker scattering of the quasiparticles in the 
BSCCO single-crystal whiskers used in our experiments. 
At higher temperatures the data are also consistent with 
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FIG. 4: Magnetic field dependence of the Josephson flux-flow 
resistance, Rff at different temperatures with fits to Eq. 
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FIG. 6: Solid triangles show temperature dependencies of 
the out of plane quasiparticle conductivity a c (a) and in- 
plane quasiparticle conductivity a a t (b). Below T c they ex- 
tracted from the JFF experiment on BSCCO long stack #4 
and above T c they represent the dc normal-state conductivi- 
ties of whiskers measured independently on samples from the 
same batch. Open circles correspond to the a c data from 
Ref.[H obtained on small mesas in zero field, open squares 
correspond to 14.4 GHz microwave data for a a t, from Ref. ITU 
obtained on BSCCO epitaxial films. Solid lines in both plots 
are just guides to the eye. Inset in the upper figure shows the 
low-temperature part of cr c (T) plotted versus T 2 . 



the results of DC measurements of a a b{T) in the nor- 
mal state carried out at the whiskers from the same 
batch^ (see the points at T > 100 K). We see that, 
in contradiction with the naive SCTMA predictions, 5 
cr a b(T)/a ab (0)-l ~ (Tf-ft) 2 , a ab (T) has strong tempera- 
ture dependence at T < ji ss 30 K so that the saturation 



to the low-temperature value is not reached even at 4 K. 

In contrast, the temperature dependence of a c ex- 
tracted from the flux-flow experiment (Fig. is con- 
sistent with cr c (T) measured on small mesas in zero mag- 
netic fieldfii with <7 C (0) being close to the universal value 
cr c (0) w 2 (kfi cm)" 1 predicted in Ref. [U The data 
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FIG. 7: Temperature dependence of the parameter V 
(cr a b<Jc)/'y 2 . Solid line is a guide to eyes. 



also well extrapolate to the points of dc measurements of 
<7 C (T) of the same mesas in the normal state. The tem- 
perature dependence of c c (T) is in agreement with the 
theoretical prediction Eq. (JJJ with 7, s» 23.5 K. 

Clearly, behavior of a a b(T) and cr c (T) is very differ- 
ent. This can be seen more clear in the temperature 
dependence of their ratio. In Fig. [7] we plot temperature 
dependence of the parameter T = {<J a b / & c) ll 2 (the ratio 
of the QP dissipation anisotropy to the superconduct- 
ing anisotropy). The important point is that both com- 
ponents of the quasiparticle conductivity were extracted 
from the same experiment. As was mentioned above the 
in-plane contribution to the flux-flow resistivity becomes 
considerable when r>l. Figure [7| shows that this con- 
dition is indeed valid in a wide temperature range below 
T c . Note that the value of V approaches 1 when T ap- 
proaches T c . The value of T also decreases below 20 K. 
The interesting issue is the low-temperature limit of T. 
Without Fermi-liquid and anisotropic scattering correc- 
tions to the universal electrical intralayer conductivity, T 
should approach unity at low temperatures, T <C ji, but 
obtained T is still well above unity at 4 K. 



B. Nonlinear Josephson flux-flow regime. 

A non-linear Josephson flux-flow occurs at high veloc- 
ities of Josephson vortex lattice, especially at velocities 
approaching the minimum velocity of the electromagnetic 
wave (Swihart velocity), c s = cs/(2A a (,y / e). The main 
source of the nonlinearity is the pumping the energy from 
a dc source into the travelling electromagnetic wave, gen- 
erated by the moving lattice. Due to the deformation of 
the the moving lattice induced by interaction with the 
boundaries^ the I-V characteristics in this regime can be 
calculated only numerically. Below we present the main 
steps of a calculation (see Ref. |4fJ for details) as well as a 
comparison with the experiment. We will consider here 
only the main flux-flow branch corresponding to the case 
of uniformly sliding triangular lattice. The termination 
point of this branch is related to instability of the moving 
triangular vortex lattice42^i The high-voltage branches 



in the JFF regime are apparently related with more com- 
plicated behavior and beyond the scope of this paper. We 
emphasize that at high velocity, when the washboard fre- 
quency exceeds the inverse quasiparticle relaxation time, 
one has to take into account the frequency dependence of 
the quasiparticle conductivity 17 which leads to the renor- 
malization of the plasma frequency and Swihart velocity. 

At high fields in the resistive state the interlayer phase 
differences depend approximately linearly on coordinate 
and time 



t (t, a:) « wst + k H x + <j> n 



(4) 



where u>e is the Josephson frequency and ku is magnetic 
wave vector. In the following we will use reduced param- 
eters: cue —> loe/^p, kji — + kjjjs (see Table 1). The 
most important degrees of freedom in this state are the 
phase shifts </> n , which describe the structure of the mov- 
ing Josephson vortex lattice. In particular, for the static 
triangular lattice <p n = nn. The lattice structure experi- 
ences a nontrivial evolution with increase of velocity. The 
equations for <j) n can be derived from the coupled sine- 
Gordon equations for 6 n (t,x) by expansion with respect 
to the Josephson current and averaging out fast degrees 
of freedom^ In the case of a steady state for a stack 
consisting of N junctions, this gives the following set of 
equations 



- ^2 l m [G ( n i m ) ex P 1 

m— 1 



(t>n))\ = ij; 



(5) 



where n = 1, 2 ... N and ij = ij(kn, ^>e) — (sin# n (i, x)) 
is the reduced Josephson current, which has to be ob- 
tained as a solution of these equations. The complex 
function Q{n, m) describes phase oscillations in the m-th 
layer excited by the oscillating Josephson current in the 
n-th layer. For a finite system it consists of the bulk term 
G(n — m) plus the top and bottom reflections (multiple 
reflections can be neglected): 

g(n,m) = G(n-m)+BG(n+m)+BG(2N+2-n-m), 



where 



G(n) 



dq 



exp (iqn) 



n 2 {k,q,uj) 



2tt lu 2 — iv c LO — Q 2 (k, q, u>) ' 

k 2 (1 + iVgbUj) 

2(1 - cos q) + (1 + iv ab w) /I 2 



(0) 



u> = loe and k = kn are the frequency and the in-plane 
wave vector of the travelling electromagnetic wave gen- 
erated by moving lattice. The real and imaginary parts 
of £l 2 (k,q,uj) give the spectrum of the collective plasma 
oscillations and their dumping due to in-plane quasipar- 
ticle dissipation (in reduced units). The dissipation pa- 
rameters, v c and v a b, and reduced penetration depth I 
arc defined in Table 1. B = B(k,u>) is the amplitude of 
reflected electromagnetic wave. For the practical case of 
the boundary between the static and moving Josephson 
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lattices a detailed calculation of B(k,u>) is presented in 
Ref . HI 

In general, the quasiparticle conductivities in the def- 
initions of v c and v a b are the complex conductivities at 
the Josephson frequency. The frequency dependence is 
especially important for the in-plane conductivity. The 
maximum Josephson frequency at the termination point 
of the flux- flow branch exceeds the value 1/r at fields 
B > 2 T. Therefore the frequency dependence has to 
be taken into account. We use the Drude-like frequency 
dependence of v a b = v a b(u): 

where r is the quasiparticle relaxation time. The solution 
of Eq. ||SJ| yields the current-voltage characteristic 

j(E z ) = a c E z + J c ij(k Hl u E ), 

where fc# and ug has to be expressed via magnetic and 
electric fields (see Table 1). 

We solved Eqs. numerically and calculated the I- 
V dependencies for the first flux-flow branch. We used 
<7 C and CTah/7 4 obtained from the fit of pff(B), assumed 
X a b = 200 nm, and adjusted 7 to obtain the I-V depen- 
dencies most close to experimental ones. The results are 
shown in Fig.[S]for two values of magnetic field, B = 1 T 
and B = 2 T at T = 4.2 K. One can see a very reasonable 
fit to experiment for both field values using 7 « 500. At 
high fields (see curves at 2 T) we found the fit can be 
significantly improved by taking into account frequency 
dependence of a a b via Eq. Q. The best approximation 
here was found for l/(27rr) = 0.6 THz. Both fitting pa- 
rameters found to be quite reasonable. The value for 7 
is consistent with the typical value of J c (0) for our sam- 
ples at low temperatures, J c (0) = 1 — 2 A/cm 2 , while 
the value of 1/ (27tt) is consistent with the microwave 
data of Corson et alM* where the relaxation rate at low 
temperature was found to be ~ 1 THz. 



C. Renormalization of Swihart velocity by 
quasiparticles 

Relatively low quasiparticle relaxation rates lead to an 
important observable consequence: renormalization by 
quasiparticles of the Swihart velocity at high frequencies, 
u > 1/t. The Swihart velocity is the in-plane velocity 
of the electromagnetic wave at the maximum c-axis wave 
vector, k z = tt/s. For this mode, out-of-phase oscillations 
of c-axis current in the neighboring layers induce oscilla- 
tions of the in-plane current leading to strong coupling 
with the in-plane charge transport. As a consequence, 
the velocity of this mode, cs/(2A a &v / e), is strongly re- 
duced by the Cooper pairs in comparison with veloc- 
ity of the transverse electromagnetic wave, c/y/e. At 
small frequencies the in-plane motion of quasiparticles 
only contributes to the dissipation. However, when fre- 
quency exceeds the scattering rate, quasiparticles con- 
tribute to the inductive response in the same way as 
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FIG. 8: Comparison of the experimental I-V characteristics 
measured on long BSCCO stack at fields B = IT and B = 2T 
at 4.2 K with calculated numerically using Eq. @. For both 
cases 7 was used as a fitting parameter. The best fit was 
found at 7 = 500. For the I-V characteristic at B=2 T also 
t has been varied. The best fit was found at lot — 0.25. 
For such relaxation rate frequency dependence of a a b have 
very weak influence on the I-V shape at 1 T. In the upper 
axes the voltage scale V mia (B) = Bs 2 / (2\ aby /e) ~ 0.5 mV/T 
corresponds to the voltage at which the lattice velocity reaches 
the minimum velocity of electromagnetic wave sc/ (2\ a b\ft)- 
The parameters assumed in calculations are listed in the lower 
plot. 



the superconducting electrons do, reducing an effective 
screening length and increasing the Swihart velocity. We 
will analyze this effect quantitatively. The spectrum of 
the plasma mode, oj p — u> p (k, q) and its damping param- 
eter v — v(k, q) are determined by the equation (reduced 
iv c Lu p + fi 2 (fc, g, ujp) with Cl(k, q, u p ) 
i.e., u p (k, q) is given by the solution 



units) u^ + ivUp = 
defined by Eq. ©, 
of the equation 



-Re[n 2 (k,q,u p )} 

For the minimum frequency at fixed k corresponding to 
q = ir we obtain 

T 



l + Re 



1 + ITUp 



From this equation one can easily observe that the Swi- 
hart velocity c 8 = u p /k (in units of sc/(X a b\/e)) has two 
simple asymptotics 

1/2, at c s kr < 1 



V 1 + v a b/T/2, at c s kr > 1 
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The physical meaning of the quasiparticle renormaliza- 
tion factor r = y/l + v a b I t becomes more transparent 
if we transfer to the real units and use the two-fluid ex- 
pressions for X ab and a ab , \ 2 ab = m ab c 2 / (Aim s e 2 ), a a b = 
n n Te 2 /rriab- Then the renormalization factor reduces to 
v a b It — n n /n s and 

f J irs 2 n s e 2 /erriab, at c s kr <C 1 

c s = s 

Is/its 2 (n s + n n ) e 2 /sm a b, at c s kr > 1 

i.e., the renormalization amounts to replacement of the 
superfluid density n s by the total density n s + n n . 

The first flux-flow branch terminates due to instability 
of the triangular lattice configuration at the velocity close 
to the Swihart velocity42^i Therefore, the Swihart veloc- 
ity can be estimated from the voltage at the endpoint of 
the first flux-flow branch. As a consequence of quasipar- 
ticle renormalization, one can expect an increase of the 
maximum voltage for the first flux-flow branch when the 
Josephson frequency at this voltage exceeds 1/r. The 
estimate of the renormalization factor for our samples 
for ujpT = 0.25, A a6 = 0.2 /im, a ab = 4 • 10 4 (fi cm) -1 
gives the value r = 2. Experimentally, the Swihart veloc- 
ity can be extracted from the slope of the linear depen- 
dence of the maximum flux-flow voltage on the magnetic 
field, V ffmax (H)/H. We found (see also Ref.HJ, that 
the slope increases at magnetic fields above 1.5 T (Fig. 
|SJ. That field corresponds to the washboard frequency 
about 0.5 THz, which is very close to the quasiparticle 
relaxation rate. Therefore, the increase of the slope of 
Vf fmax(H) at higher frequencies may be interpreted as 
an increase of the Swihart velocity due to a renormaliza- 
tion of plasma frequency at u> > 1/r. Experimentally, the 
slope of Vff m ax(H) increases by factor 2.1 that is very 
close to the theoretical estimate of the renormalization 
factor. 



IV. CONCLUSIONS 

In summary, we carried out detailed experimental and 
theoretical studies of the linear and non-linear Joseph- 
son flux-flow regimes in BSCCO. Both regimes as shown 
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Table 1. Meanings, definitions and practical formulas for the reduced parameters used in the paper. In practical 
formulas f p = lo p /2tt means plasma frequency, p c and p a }, are the components of the quasiparticle resistivity 



Notation 


Meaning 
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Prar-tir-al fnrmula fRSPPnl 

_riacLicai loimuia ^_dov_^ov^j 
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magnetic wave vector 
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c-axis dissipation parameter 


4ir<r c 
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in-plane dissipation parameter 
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c 2 
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reduced London penetration depth 
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